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Refractive  index  and  wavenumber  properties  for  cyclotron  resonant 
quasilinear  diffusion  by  cold  plasma  waves 
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Wave-particle  interactions  have  a  large  effect  on  magetospheric  particles,  in  the  radiation  belts  and 
elsewhere.  Bounce-averaged  quasilinear  diffusion  coefficients  have  been  calculated  for  whistler  hiss 
and  chorus  and  electromagnetic  ion  cyclotron  waves  (EMIC),  which  are  all  believed  to  play  major 
roles.  To  perform  these  calculations  efficiently,  techniques  have  been  developed  that  use  properties 
of  the  refractive  index  of  these  modes  to  identify  ranges  of  wave-normal  angle  that  are  compatible 
with  cyclotron  resonance  in  a  given  frequency  band.  Other  cold  plasma  waves,  in  the  L-X,  L-O. 

R-X,  and  Z  modes,  can  also  resonate  with  energetic  electrons,  and  some  preliminary  calculations  of 
their  diffusion  coefficients  have  been  reported.  Here,  it  is  shown  that  the  refractive  index  of  these 
modes  allows  the  techniques  developed  for  whistler  and  EMIC  waves  to  be  applied  to  them  as  well. 
Sample  calculations  are  presented  for  Z  mode  waves,  with  a) ^  >  f le  and  tope  <  (le.  It  is  also  observed 
that  for  any  cold  plasma  mode,  the  wavenumber  is  an  increasing  function  of  frequency  for  a  fixed 
value  of  wave-normal  angle;  this  is  proved  algebraically  with  mild  approximations  and  verified 
numerically  for  a  very  wide  range  of  parameters.  This  allows  a  variant  of  the  technique  for 
efficiently  calculating  diffusion  coefficients.  [DOI:  10.1063/1.2744363] 


I.  INTRODUCTION 

Cyclotron-resonant  wave-particle  interactions  are  crucial 
to  the  dynamics  of  magnetospheric  particles,  particularly  ra¬ 
diation  belt  electrons.  Whistler  mode  hiss  w'ithin  the  plasma- 
sphere,  whistler  mode  chorus  outside  the  plasmasphere,  and 
electromagnetic  ion  cyclotron  (EMIC)  waves  at  the  plasma- 
pause  or  within  density  plumes  are  all  believed  to  play  major 
roles.  Their  quantitative  effects  have  recently  been  calculated 
in  detail, 1-0  using  bounce-averaged  local  quasilinear 
theory.7  8  In  the  low-density  region  outside  the  plasmasphere, 
electrons  are  also  capable  of  resonant  interactions  with  the 
cold  plasma  L-X,  L-O,  R-X,  and  Z  modes.9,10  Preliminary 
calculations  of  diffusion  coefficients  for  these  modes  have 
also  recently  been  presented.5,11 

The  full  bounce-averaged  quasilinear  calculations  are 
fairly  demanding,  since  even  local  diffusion  coefficients  in¬ 
volve  integration  over  a  continuum  of  waves,  requiring  the 
resonant  frequency  for  each  wave-normal  angle  or  vice 
versa,  compounded  by  the  sum  over  resonant  harmonic  num¬ 
ber  n  (in  principle,  from  to  °°).  When  bounce  averaging  is 
performed,  it  becomes  impractical  to  repeat  the  calculations 
for  many  different  values  of  particle  energy  and  equatorial 
pitch  angle,  background  density,  and  wave  parameters.  To 
circumvent  this,  simplified  forms  of  the  local  diffusion  coef¬ 
ficients  have  been  derived  assuming  parallel  propagating 
waves,  *  and  bounce  averaged  and  found  to  give  good  esti¬ 
mates  of  the  exact  results  for  chorus  waves  in  several  test 
cases.'  A  different  approach  is  to  identify  in  advance,  and 
eliminate,  integration  ranges  that  are  incompatible  with  cy¬ 
clotron  resonance  for  the  given  wave  population.  This 
greatly  reduces  the  amount  of  computation  required  without 
changing  the  result,  provides  a  rational  criterion  for  limiting 
the  range  of  /7,  and  provides  insight  into  which  waves  con¬ 
tribute  to  the  diffusion  (which  suggests  an  approach  for  pa¬ 


rameter  sensitivity  studies).  Such  techniques  have  previously 
been  developed  for  whistler  mode  hiss,1  EMIC,2  and  chorus,3 
based  on  detailed  analysis  of  the  refractive  index  of  these 
wave  modes.  Here,  it  is  shown  how  they  may  be  easily  ap¬ 
plied  to  L-X,  L-O,  R-X,  and  Z  mode  electromagnetic  waves 
as  well. 

Also  discussed  is  a  closely  related  approach  based  on  the 
behavior  of  the  wavenumber,  k ,  rather  than  refractive  index, 
kef  a >,  for  cold  plasma  waves.  It  is  shown  in  general  that  k  is 
an  increasing  function  of  frequency  o),  at  a  fixed  wave- 
normal  angle  0.  This  basic  observation  about  k(io)  does  not 
seem  to  have  received  as  much  attention  as  k(0)  has.  Here 
this  property  is  established  theoretically  as  far  as  possible, 
and  numerically  in  general,  to  justify  its  use  in  calculating 
diffusion  coefficients. 

II.  WAVE  MODES 

All  of  the  cold  plasma  wave  modes  are  given  by  the 
standard  expression  for  the  refractive  index,14 

2 

^±=~  =  {(RL  -  PS) sin2  6+  IPS 
k±c~ 

±  a \(RL-PS)2  sin4  0+4P7D2  cos2  0)/2PRL , 

(1) 

where  the  quantities  (R%L,PyS,D)  depend  on  the  wave,  cy¬ 
clotron,  and  plasma  frequencies  (w,  fl,  and  cop,  respectively), 
and  (7  is  defined  to  be  +1  for  PDX)  and  -1  for  PD<().  The 
to  dependence  for  a  pure  hydrogen  plasma  is  shown  in  Fig.  1 
for  u)2pe/£l*= 0.3  and  u)2peliVe  =  2.2  (representative  of  back¬ 
ground  plasma  with  mildly  low  and  mildly  high  density,  re¬ 
spectively),  for  several  values  of  0 .  Also  shown  are  the  lim¬ 
iting  cases  >I'_(0=O)=1//?,  ^+(0=0)=  1  !  Ly  and  ¥+(0 
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FIG.  1  The  functions  V ±= «r / k^c2  vs  a >/Cle  for  a  pure  hydrogen  plasma,  with  0=  15°  ,30°  ,45°  ,60°  ,75°.  The  superluminous  L-X,  L-O,  and  R-X  modes,  and 
the  whistler  and  Z  modes  are  shown  with  solid  curves  for  the  branches  and  dashed  curves  for  yi/_.  The  thick  curves  with  J  are  1  !L,  MP,  and  MR  (left 
to  right,  respectively);  the  thick  curve  with  ^<1  is  also  MR.  The  insets  show  vs  ta/1^  for  the  same  parameters;  the  thick  curve  is  ML,  with  maximum 
value  at  ta=0.  For  each  wave  mode  in  the  top  panel,  the  arrows  show  the  direction  of  increasing  0. 


=  90°)=  1  IP,  as  well  as  the  frequencies  o)L ,  a)P=(ope,  and  coR 
at  which  L,  P ,  and  R  are  zero,  respectively.  The  branches 
above  the  line  (o2/k2c2=  1  are  superluminous,  with  phase  ve¬ 
locity  cj/k>c.  Finally,  the  insets  show  VP+  for  0<o><f2,, 
which  are  EMIC  waves. 

With  the  heavy-ion  or  magneto-ionic  approximation 
melmx- 0  in  the  definitions  of  R ,  L,  and  P,  Eq.  (1)  can  be 
written  in  the  form 


't'±  =  [l-l/(/t<D2A±)]-\ 

, _  (2) 

^  ^  sin2  0  j  A2  sin4  0  cos2  0 

±_  ~  2 Mr  -  1)  ±  V  4 (A<b2  -  l)2  +  w2  ’ 


where  a>=a >/ile  and  A—Q2/ (o2.  Also,  o)RL 
=  ^e[(l+4M)1/2±l]/2  and  cjp=i}r/ A]l .  Note  that  (oL<wP 
<a>R,  and  that  coL>il,  as  long  as  <o2pcl£l2e>mef  m,.  Positive 


'P  requires  either  Acj2 A>1  (superluminous  modes)  or  A 
<0  (subluminous  modes). 

The  three  superluminous  modes  are  identified10  as  the 
L-X,  with  (i)L<u)<u)P  and  described  by  ^P+;  the  L-O,  with 
u)>u)P  and  described  by  'P*;  and  the  R-X,  with  (o>  cjr  and 
described  by  y¥_.  Also  shown  are  the  subluminous  whistler 
and  Z  modes,  given  by  ^P_.  The  whistler  mode  has  co 
<min[ft<,,(u/>],  while  the  Z  mode  propagates  for  oj  between 
<x)P  and  the  upper  hybrid  frequency  o>UH.  defined  by  5=0. 
This  lies  below  <oR  since,  in  the  magneto-ionic  approxima¬ 
tion,  <oim=(le(]  +  \/A)m  so  that  a>^-6j2H=^A>0. 

Both  whistler  and  Z-mode  propagation  can  be  subject  to 
a  resonance  cone,  given  by  tan20RC=-P/Ss=(co2-a)2)(fl2 
-co2)/(o^jH-(o2)(o2,  assuming  (0>(oLHs*({lifle)I/\  Propaga¬ 
tion  then  requires  0<  0RQ.  Z-modc  waves  will  have  a  reso¬ 
nance  cone  when  <u>fl,,  requiring  0>  ^  for  propagation; 
this  is  necessarily  the  case  when 
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FIG.  2.  The  functions  'P  (solid  curves)  and  V  (dashed  curves)  with  o>^/flj=0.3,  wave-normal  angle  0  from  0  to  90°,  and  pitch  angle  a  from  0  to  90°,  for 
«*-l,  -2.  and  +1.  Cyclotron  resonance  occurs  w'here 


III.  DIFFUSION  COEFFICIENTS  AND  RESONANCES 

To  implement  the  quasilinear  formulation,  the  wave 
population  is  usually  modeled  by  the  product  of  Gaussian 
distributions  of  frequency  and  wave-normal  angle  (0  or  jc 
=  tan  0),  with  sharp  cutoffs  for  each.7,8  The  local  diffusion 
coefficients  (for  pitch  angle,  momentum,  and  cross  diffusion) 
have  the  form  of  sums  of  integrals  over  0  from  a  minimum  to 
a  maximum  value. 


keep  £le  and  oope  fixed,  it  may  be  useful  to  redraw1^  the  CM  A 
diagram  delineating  the  various  wave  modes,  with  axes 
w2r/H2  and  <j)t (le  rather  than  the  conventional  quantities 
co^/co2  and  £lef  co). 

The  cyclotron  resonance  condition  can  be  written  as 


Ire2  (co  -  co„)2  c1 


cos2  a  cos2  6  =  V 


(5) 


D  = 


'Dn(Q,o))dQ. 


(3) 


Here  go  is  the  resonant  frequency  corresponding  to  n  and  0 , 
given  by 

6>  -  =  o)fn  ion  =  sniljy,  (4) 

where  12c  is  the  nonrelativistic,  unsigned  cyclotron  frequency 
of  the  diffusing  test  particle  and  s=q/\q\.  Vn  is  taken  to  be 
zero  unless  go  lies  between  lower  and  upper  cutoffs  <oLC  to 
o>lJC.  Within  the  cutoffs,  the  distributions  may  be  param¬ 
etrized  by  peak  and  half-width  values/  e.g.,  (0m,  0>v.;o>m,  <Sco). 
The  permissible  frequency  range  depends  on  the  wave  mode 
being  studied,  but  the  range  actually  used  constitutes  an  em¬ 
pirical  model  of  the  waves  and  may  be  much  narrower. 
(Since  the  local  diffusion  coefficients  involve  ranges  of  co  but 


This  includes  resonances  with  both  signs  of  which  is 
justified  if  the  upgoing  and  downgoing  wave  populations  arc 
assumed  equal,  or  if  both  halves  of  the  particle’s  bounce 
motion  are  considered. 

Figure  2  shows  V(o>)  and  ^(w)  for  co2pe/iVf=03,  for  a 
wide  range  of  6  and  a.  The  'P  curves  are  the  same  in  each 
panel  [and  the  same  as  in  Fig.  1(a)],  while  F,  which  de¬ 
creases  as  6  increases,  is  shown  for  different  values  of  en¬ 
ergy  and  n.  Generally  speaking,  with  n  =  - 1  resonances  arc 
readily  found  for  the  L-X  and  L-O  modes  (as  well  as  the 
whistler  and  Z  modes)  but  are  marginal  for  the  R-X  mode,  in 
agreement  with  Ref.  10.  However,  with  n  =  - 2.  resonances 
are  also  plentiful  for  the  R-X  mode  over  a  broad  range  of 
parameters,  as  are  resonances  with  n>0  for  the  whistler  and 
Z  modes. 

To  recap  the  calculation  method  of  Refs.  1-3.  the  idea  is 
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to  identify  the  ranges  of  0  within  0min  to  0Iliax  for  which  co 
lies  outside  (oLC  to  <0^,  and  exclude  them  before  performing 
the  integrations.  To  this  end,  each  side  of  Eq.  (5)  is  consid¬ 
ered  geometrically  as  a  function  of  a)  at  fixed  0.  The  mini¬ 
mum  and  maximum  values  of  V(o>)  over  (UjC  to  couc  are 
readily  determined.  For  example,  with  cun<0,  is  mono- 
tonically  increasing  so  Vmin=V((oLC)  and  Vmax= V(a>uc).  If 
the  minimum  and  maximum  values  of  'P  can  be  character¬ 
ized  as  well,  then  0  ranges  are  disqualified  according  to  the 
simple  criteria 

V  >  'P  or  V  <  'P  (6) 

Ym\n  x  max  Ul  *max  Tmin*  \u' 

These  conditions  are  easily  interpreted  graphically  by  draw¬ 
ing  a  rectangle  defined  by  (o>lc^a>^o>uc,  'Pmin<'vP 
The  P(tu)  and  ^(<0)  curves  can  only  intersect  for  6 
such  that  V  passes  through  this  rectangle;  other  ranges  of  6 
need  not  be  considered.  An  extreme  instance  of  this  is  that 
superluminous  modes  have  ^^>1  while  Vmax<l  for  o)n 
<0,  so  there  are  no  superluminous  resonances  with 
(i.e.,  with  n^0  for  electrons).  Algebraically,  conditions  (6) 
are  easily  solved  quadratic  inequalities1’^  for  sin2 0.  Exclud¬ 
ing  the  corresponding  6  ranges  greatly  reduces  the  amount  of 
work  required  to  evaluate  the  diffusion  coefficients,  particu¬ 
larly  since  for  \n\  sufficiently  large,  no  6  values  remain  at  all. 

The  success  of  the  method  depends  on  being  able  to  find 
(or  bound)  the  minimum  and  maximum  values  of  'P(aj)  for 
fixed  0.  This  is  especially  simple  if  ¥( co )  behaves  monotoni- 
cally,  as  Fig.  1  suggests  often  occurs.  As  noted  above,  this 
has  already  been  studied  in  detail  for  chorus  mode  and  EMIC 
waves.  In  the  following  subsections,  the  behavior  of  'P(o))  is 
investigated  for  the  superluminous  and  Z  modes. 

A.  Superluminous  modes 

From  Eq.  (2),  the  sign  of  the  slope  'P'  =  <7vP/<? o>  for  the 
superluminous  modes  is  given  by 

^--(orA^^Aj.  (7) 

Then  it  is  straightforward  to  show  that  ^<0,  so  that  'P' 
<  0  for  both  the  L-X  and  L-O  modes.  For  the  R-X  mode,  it 
can  be  shown  that  A^.>0  for  cj><or  (or  even  for  the  less 
demanding  condition  (o>coP ).  It  is  also  easily  checked  that 
A_>0  at  ojr ,  and  therefore  at  co>gjr,  so  that  'P'<0  holds 
for  the  R-X  mode  as  well.  Thus  conditions  (6)  may  be  ex¬ 
ploited  for  superluminous  waves  just  as  it  was  for  EMIC 
waves." 

B.  Z  mode 

Figure  3  shows  'P  for  the  whistler  and  Z  modes  with 
iope/ilj= 0.1  and  0 =  1  0  ,2°  ,3°.  It  is  evident  that  the  Z  mode 
does  not  strictly  obey  ^'(cu)>0  for  a)<Cle/ 2,  where 
\/R(co)  is  increasing,  if  0  is  small  enough  that  'P—l/ft. 
Such  Z  mode  waves  require  o)pf<u)<CltJ 2,  or  a)pe/£le 
<0.5.  It  is  conjectured  that  this  behavior  is  controlled  by  the 
sign  of  yV"  evaluated  at  fl^/2,  just  as  'P"  evaluated  at  <o=0 
was  seen3  to  control  the  behavior  of  the  whistler  mode  for 
very  small  oj  and  large  0.  (Perhaps  fortuitously,  the  numeri¬ 
cal  example  of  Ref.  5  takes  ojpe/£le=Q.5  and  g>lc/1\=0.55). 


Wp 


FIG.  3.  The  function  Nk  ._= <u2 / Ire2  vs  w/ile  for  the  whistler  and  Z  modes, 
with  fP/ a/5  =0.1  and  0=1°, 2°, 3°.  Also  shown  is  the  value  for  0=0, 
namely  1  /R. 


1.  Estimates  of  *I'max  and  ^Pm,n 

To  characterize  the  behavior  of  ^(tu)  in  the  problematic 
cases,  we  turn  to  a  factorization  found  for  the  whistler  mode, 
'P  =  ‘VP1'P2>  where 


D  <0 


V',= 


nj  (RL  -  PS)sin2  0+2  PS 


coL  [ 


2  PD 


(RL-  PS)2  sin4  0 


4P2D2 


+  COS2  0  I 


(8) 


For  (o2>flrfli,  S',  is  an  increasing  function  of  to  while  S',  is 
a  decreasing  function  of  o>;  a  proof  was  outlined  in  Ref.  16. 
(Walker,17  in  a  related  context,  refers  to  “the  horrendous  al¬ 
gebra  involved/')  A  careful  review  of  that  proof  shows  that  it 
holds  for  the  Z  mode  as  well.  Thus  the  bounds 


may  be  exploited  for  Z  mode  waves  just  as  it  recently  was 
for  whistler  waves.3  Actually,  if  ^LC<^?/-<a>uo  treating 
the  two  subintervals  separately  gives  the  tighter  bounds 

min[^Pj (t^Lc)^/2(^f/^)»^/(tA,uc)]  ^  ^P(<o) 

<  maxt'P^tuuc)^^^), ^(11^/2)].  (10) 


Thus  one  may  safely  skip  0  ranges  for  which  both  Vrmin 
>^]((ouc)‘vP2(a>|  c)  anc*  V/min> ^(11^/2)  hold,  as  well  as  0 
ranges  which  satisfy  both  Vmax<vP1(6uI ^^(a^c)  and 

Yet  another  estimate  of  *Pmin  and  xPn)ax  for  the  Z  mode 
comes  from  the  observation  that  xJ/(<o)  increases  with  in¬ 
creasing  0.  To  prove  this,  one  can  show  in  the  magneto-ionic 
approximation  that  XP1>0  for  u)<coR  and  that  rAp2/00  has 
the  same  sign  as  L  does,  namely  positive  for  the  Z  mode. 
Thus  'P  is  bounded  by  min[vP(w,  (9min)]  and 
max[xP(u>,  0max)],  which  may  be  found  numerically  and  ex¬ 
ploited  in  Eq.  (6). 
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FIG.  4.  Diffusion  coefficients  for  10  pT  Z-mode  waves  within  10°  of  the  equator.  The  top  row  shows  equatorial  and  bounce-averaged  pitch  angle  (solid 
curves),  momentum  (dashed),  and  cross  (dotted)  diffusion  coefficients  for  £le=0.5  at  the  equator.  The  local  values  generally  lie  above  the  corresponding 
bounce-averaged  quantities.  The  bottom  row  shows  bounce-averaged  diffusion  coefficients  for  u)pe/{l,=4  at  the  equator. 


2.  Numerical  results 

Figure  4  shows  diffusion  coefficients  for  electrons  inter¬ 
acting  with  Z-mode  waves  at  L=4.5.  The  same  wave  model 
parameters  are  used  as  for  the  local  equatorial  calculations  of 
Ref.  5,  namely  a>/v/flc  =  0.5,  with  frequency  parameters 
( ,  <$6> ,  o)LC ,  <ouc)  =  (0.60 , 0.05 , 0.55 , 0.65)12,  and  wave- 
normal  angle  parameters  (0m,  0W,  0min,  0mJ 
=(0,30°  ,0,45°),  as  well  as  wave  amplitude  10  pT.  The  top 
row  of  Fig.  4  shows  both  local  equatorial  and  bounce- 
averaged  diffusion  coefficients;  the  latter  took  Z-mode  waves 
to  be  present  only  for  latitude  within  10°  of  the  equator  (with 
constant  u)pe  and  dipolar  variation  of  11,).  The  local  values 
replicate  the  values  of  Ref.  5  (when  appropriately  scaled  to 
1 00  pT),  while  the  bounce-averaged  results  are  significantly 
smaller.  For  0.1  MeV,  is  less  than  Dpp  but  for  higher 
energy  the  reverse  is  generally  true.  Also,  the  waves  seem 
ineffective  for  intermediate  pitch  angles,5  with  implications 
for  the  resulting  particle  distribution.  However,  such 
conclusions  depend  on  the  parameter  choices,  as 
shown  in  the  bottom  row  of  Fig.  4.  Here  (°peltoe= 4, 
with  (ojm ,  Sco ,  o>uc ,  c)  =  (4.02 ,1.0,4.01. 4.025)11,  and 
0mi,n  0max)  =  (5O°  ,50°  ,30°  ,80°).  The  frequency 
cutoffs  span  about  140  Hz,  centered  near  38  kHz.  For  larger 
<op,/12„  the  restrictions  on  frequency  and  wave-normal  angle 


start  to  become  prohibitive.  Nevertheless,  the  diffusion  coef¬ 
ficients  shown  are  substantial  over  a  large  range  of  pitch 
angle,  and  the  momentum  rates  can  be  larger  than  the  pitch 
angle  rates,  a  feature  expected  at  low  pitch  angle  for  waves 
with  o>>12,.8  For  the  amplitude  assumed,  the  rates  are  com¬ 
parable  to  values  computed  for  other  waves, ig  though  local¬ 
ization  of  the  waves  in  longitude  has  not  been  accounted  for. 
Apparently  Z-mode  waves,  if  present  under  favorable  condi¬ 
tions,  are  capable  of  resonating  effectively  with  radiation  belt 
electrons. 

IV.  VARIANT  OF  THE  METHOD  BASED 
ON  WAVENUMBER 

The  resonance  condition  may  also  be  written  as 

kc  =  |cu  -  coJ(cA;)sec  a  sec  0  =  12,  W,  (11) 

again  keeping  both  signs  of  k\V\\.  (At  a  fixed  location,  12,  is 
constant  and  just  serves  to  normalize  W.)  Figure  5,  analogous 
to  Fig.  1,  suggests  that  k ,  or  A2c2=o>2/' I',  is  an  increasing 
function  of  (o  at  fixed  0  for  all  of  the  cold  plasma  wave 
modes.  This  is  supported  below.  [This  monotonic  behavior 
was  shown  in  a  similar,  comprehensive  set  of  plots  by 
Walker18  (Figs.  3-18  to  3-20)  but  not  remarked  upon,  al¬ 
though  “contemplation”  of  these  and  related  figures  was  rec- 
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a  l  up  aH 


FIG.  5.  The  functions  k^c2lil^  vs  a>/(le  as  in  Fig.  1(a).  The  diagonal  line  separates  the  superluininous  from  the  subluminous  modes. 


ommcndcd.]  This  behavior  immediately  suggests  conditions 
analogous  to  those  in  Eq.  (6)  to  rule  out  resonances  in  a 
given  frequency  range, 

^min  >  (kc)mdLX/Cle  or  VVmax  <  (&c)mir/Qr,  (12) 

where  (kc)min  is  evaluated  at  coLC  and  (kc)max  is  evaluated  at 
£uuc.  As  before,  this  leads  to  simple  quadratic  equations  for 
sin2  0.  The  following  subsections  verify  the  property  (k2)' 
>0  for  the  various  wave  modes. 

A.  Superluminous  and  EMIC  waves 

For  the  superluminous  modes  in  the  magneto-ionic  ap¬ 
proximation,  co2/k2c2  has  been  shown  to  be  a  decreasing 
function  of  co  so  obviously  \lk2c2  is  a  decreasing  function  of 
co  and  k2c2  is  increasing.  The  same  holds  for  EMIC  waves, 
even  in  multispecies  plasmas,  in  the  high-density 
approximation.2 

B.  Z  mode 

For  the  Z  mode,  the  sign  of  the  slope  of  k2  is  given  by 

(klY  ~2Ao>2A2  +  A:  (13) 

and,  as  noted  for  the  R-X  mode  in  Sec.  Ill  A,  A^  >0  when 
A  co2  -  1  >0  (which  is  guaranteed  by  co>coP).  Thus  (k2)'  > 0 
for  the  Z  mode. 


C.  Whistler  mode 

To  analyze  {k2)'  for  the  whistler  mode,  we  resort  to  the 
factorization  discussed  in  Sec.  Ill  B,  or  k2c 2 

=  (<i>2/'I',)/'I'2.  Since  ^2  *s  a  decreasing  function  of  co  for 
co2>nililei  to  show  that  k2  is  increasing  it  is  sufficient  to 
show  that  cu2/'!7!  is  increasing.  This  expression  is  just  coR/D , 
which  is  indeed  found  to  have  positive  slope  given  the  whis¬ 
tler  mode  properties  and  co<coP. 

D.  Numerical  tests 

Cases  not  covered  above  include  R  mode  (whistler) 
waves  with  co2<Clt,klj  and  EMIC  waves  without  the  high- 
density  approximation.  It  is  worth  recalling  the  remarks  of 
Stix 14  that  “it  is  not  difficult  for  a  problem  in  plasma  waves 
to  become  completely  unwieldy.  Even  the  simple  case  of 
waves  in  a  cold  drift-free  homogeneous  plasma  is  on  the 
borderline  of  tractability,”  and  although  “the  difficult  prob¬ 
lems  in  plasma  waves  in  a  homogeneous  medium  become, 
for  the  most  part,  just  difficult  problems  in  algebra,...  the 
algebra  is  not  trivial.”  Even  the  limiting  case  0=0  has  pit- 
falls:  (klc2)' =(co2L)'  can  in  fact  be  negative  for  il,<co 
<cuz,  but  no  such  waves  propagate  because  L  itself  is  nega¬ 
tive  in  that  range. 

To  test  the  observation  (k2Y  >0  in  full  generality,  exten¬ 
sive  numerical  tests  were  performed  as  follows.  Both  k\  and 
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FIG.  6.  Illustrations  of  two  versions  of  the  no-resonanee  criteria,  for  whistler  waves  and  600  keV  electrons;  other  parameters  are  given  in  the  text.  Resonances 
occur  where  the  curves  intersect,  (a)  Because  the  curve  of  V[u))  (dashed)  does  not  cross  the  rectangle  determined  by  the  minimum  and  maximum  values  of 
'F  (solid  curve),  conditions  (6)  are  met,  indicating  no  resonance  within  the  frequency  cutoffs,  (b)  For  the  same  parameters  as  (a),  the  curve  of  W(o>)  (dashed) 
does  cross  the  rectangle  determined  by  the  minimum  and  maximum  values  of  kc/tl,  (solid  curve),  so  conditions  (12)  do  not  rule  out  resonances  within  the 
cutoffs,  (c)  V  and  'F(w)  for  a  different  set  of  parameters.  Here,  since  the  V  curve  crosses  the  'R  rectangle,  conditions  (6)  do  not  rule  out  resonances  within  the 
cutoffs,  though  none  are  present,  (d)  For  the  same  parameters  as  (c).  W  avoids  the  kdde  rectangle,  so  conditions  (12)  successfully  identify  the  absence  of 
resonances  within  the  frequency  cutoffs. 


(£+)'  were  calculated  for  1000  values  of  A  =  £l2/ a)2pe  logarith¬ 
mically  spaced  between  0.01  and  108,  200  uniformly  spaced 
values  O^sin2  0^  1,  and  5000  uniformly  spaced  frequencies 
0<o)<2(or.  This  wide  range  of  A ,  corresponding  to  10-4 
**  <opf/{le^ \0,  includes  the  extremely  low-density  range 
u)pe<ili;  the  frequency  range  o)<H,  is  also  resolved. 
(Evaluating  these  109  cases  took  only  a  few  minutes  on  a 
moderately  fast  computer,  which  is  one  way  to  attack  the 
tractability  problem;  algebraic  manipulation  software  can 
also  be  useful.)  A  multispecies  plasma  (85%  H+,  5%  He+, 
10%  0+)  was  also  tested,  with  1000  uniformly  spaced  fre¬ 
quencies  0<co<2col.  In  every  case  with  k2>  0,  the  relation 
{k2y  >0  was  found  to  hold.  Thus  conditions  (12)  may  be 
evaluated  and  used  to  limit  the  0  ranges  of  the  integrals  for 
diffusion  coefficients. 

V.  DISCUSSION 

As  stated  above,  the  resonance  condition  may  be  written 
in  a  variety  of  ways,  and  conditions  analogous  to  Eq.  (6)  or 


Eq.  (12)  derived.  Since  all  such  methods  only  eliminate  0 
ranges  that  are  known  not  to  contain  resonances,  the  com¬ 
puted  diffusion  coefficients  should  be  unaffected,  apart  from 
considerations  of  efficiency  and  its  practical  converse,  accu¬ 
racy.  However,  variants  of  the  method  are  not  all  equally 
effective. 

Figure  6  illustrates  conditions  (6)  and  (12)  for  a  600  keV 
electron  interacting  of  whistler  waves  with  a£/ft*=0.3.  The 
plots  are  drawn  for  0=30°.  In  Figs.  6(a)  and  6(b),  the  elec¬ 
tron  has  a=25°,  the  waves  are  taken  to  have  o>,  c/ft,=0. 1 , 
a>uc/flt,=0.4,  and  resonances  with  «=1  are  considered.  In 
Fig.  6(a),  the  minimum  and  maximum  values  of  'F  for  these 
frequency  cutoffs  determine  the  shaded  rectangle.  Conditions 
(6),  equivalent  to  V  avoiding  the  rectangle,  successfully  iden¬ 
tify  the  absence  of  the  /i=  1  resonance  for  this  0.  Figure  6(b) 
shows  the  same  physical  situation,  but  conditions  (12)  for 
W{co)  to  avoid  the  rectangle  determined  by  the  minimum  and 
maximum  values  of  kc  are  not  met,  so  the  method  fails  to 
eliminate  this  nonresonant  0  value.  Figures  6(c)  and  6(d)  are 
similar,  but  with  <*=45°,  <uLC/Or=0.05,  a>lIC/n,.=0.2.  and 
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n=-\.  Here,  as  shown  in  Fig.  6(c),  the  conditions  in  (6)  are 
not  met  (the  curve  of  V  passes  through  the  Nk  rectangle),  so 
this  0  value  is  unnecessarily  kept,  while  Fig.  6(d)  shows  that 
conditions  (12)  are  successful  in  identifying  this  0  as  non- 
rcsonant  (the  curve  of  W  avoids  the  kc  rectangle). 

The  approach  should  be  most  effective  when  the  reso¬ 
nance  condition  is  written  in  a  form  f(a),0)=g(<o,0)  with 
(say)  /  as  weakly  dependent  on  o>  as  possible,  since  when  the 
curve  of  /  is  just  a  horizontal  line,  conditions  relating  the 
minimum  and  maximum  values  of  /  and  g  become  exactly 
equivalent  to  the  existence  of  resonances.  [This  is  the  case 
for  V(cj)  with  n=0,  which,  for  the  whistler  mode,  often  turns 
out  to  be  particularly  troublesome.] 

Of  course,  it  is  also  necessary  to  be  able  to  locate  the 
maximum  and  minimum  of  the  two  functions  within  coLC  to 
^uc»  w^h  0  general,  so  that  the  conditions  can  be  evaluated. 
So  far  this  has  been  done  for  the  resonance  condition  in  the 
forms  of  Eqs.  (5)  and  (11),  where  in  each  case  both  sides  of 
the  equation  are  either  monotonic  with  respect  to  at  or  else 
the  product  of  two  such  factors  (‘vJ/> j  and  Implementing 
both  versions  in  a  single  computer  code  seems  feasible  and 
worthwhile,  and  probably  also  sufficient  for  any  reasonable 
program  of  evaluating  quasil inear  diffusion  coefficients. 

Lyons7  gave  the  general  relation  between  pitch  angle, 
energy,  and  cross  diffusion  coefficients:  roughly,  energy  and 
cross  diffusion  become  competitive  with  pitch  angle  diffu¬ 
sion  when  the  resonant  frequencies  are  not  small  compared 
to  For  whistler  waves,  the  resonant  frequency  for  a  given 
particle  tends  to  increase  with  decreasing  uopc/£le ,  even 
though  this  also  decreases  the  maximum  frequency  at  which 
whistlers  can  propagate.  ’  Superluminous  and  Z-mode  waves 
are  compatible  with  (o  as  large  as  or  larger  than  £le,  which 
again  favors  energy  diffusion,  even  to  the  extent  of  exceed¬ 
ing  pitch  angle  diffusion.11  Decreasing  a)pe/ilri  especially 


below  1,  also  has  the  direct  effect  of  enlarging  the  frequency 
range  of  propagation  for  the  Z  mode,  and  Summers  cl  al.n) 
found  this  parameter  to  be  important  for  superluminous 
waves  as  well.  The  potential  importance  of  these  modes  for 
pitch  angle  and  energy  diffusion  of  radiation  belt  elections 
calls  for  more  detailed  evaluation  of  the  diffusion  rates, 
which  awaits  more  comprehensive  wave  measurements. 
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